Generating Unexpected Spin Echoes in Dipolar Solids with vr Pulses 
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NMR spin eclio measurements of ^''C in Ceo, in Y2O3, and ^^Si in silicon are shiown to defy conven- 
tional expectations when more than one tt pulse is used. Multiple vr-pulse echo trains may either freeze out or 
accelerate the decay of the signal, depending on the 7r-pulse phase. Average Hamiltonian theory, combined with 
exact quantum calculations, reveals an intrinsic cause for these coherent phenomena: the dipolar coupling has a 
many-body effect during any real, finite pulse. 
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PACS numbers: 03.65.Yz, 03.67.Lx, 76.20.+q, 76.60.Lz 

Nuclear magnetic resonance experiments rest upon a solid 
theoretical foundation 

mm 

. Advanced pulse sequences 
have been successfully applied in all subfields of magnetic 
resonance, in atomic physics, and even in the emerging field 
of quantum information processing. Thus, we were surprised 
to discover that simple experiments on doped silicon appeared 
to be inconsistent with conventional NMR theory |5]. For ex- 
ample, coherent signals may be observed well beyond the T2 
that is measured in two-pulse spin echo m\ experiments, pro- 
vided that more than one tt pulse is used [sl 0,11] • 

In this Letter, we report the same surprising phenomena 
in Buckminsterfullerene (Ceo) ™d Yttria (Y2O3), two solids 
linked to silicon through the form of the homonuclear dipolar 
coupling Ol . We also show that multiple 7r-pulse echo trains 
may either freeze out or accelerate the expected dipolar de- 
cay of the NMR signal, depending upon the phases used for 
the TT pulses. Average Hamiltonian theory i9|], combined with 
exact quantum calculations, is used to show that this pulse 
phase sensitivity has an intrinsic origin, arising from the sur- 
prisingly non-negligible effects of the dipolar coupling during 
strong but finite pulses. 

The decay of signals produced by the single-7r-pulse Hahn 
echo sequence ||6[] (HE: OOjc-r-lSOy-r-ECHO) as t is in- 
creased is a standard measure of T2 [1]. In both Ceo [Fig.[lja)] 
and Y2O3 [Fig. [Tfb)] powders, the multiple-7r-pulse Carr- 
Purcell-Meiboom-Gill sequence lO (CPMG: 90x-{r-180F- 
T-ECHO}"^P°''*) generates echoes well beyond T2. Moreover, 
CPMG in these samples show both the long tail at short r [Fig. 
[TJa) andlHb)] and the even-odd effect at long r [Fig.lTJc) and 
[Tfd)], as previously reported in silicon ijHQ]. 

At room temperature. Ceo molecules form an fee lattice, 
and each Ceo undergoes rapid isotropic rotation about its lat- 



tice point 01 in . The time-averaged spin Hamiltonian is ob- 
tained by placing ^'^C [spin I = 1/2, with 1.11% natural 
abundance (n.a.)] spins on a 60-fold degenerate fee lattice. 
Multiple occupancy of an fee lattice site is not a problem in 
this model, since the intra-Ceo dipolar couplings are averaged 
to zero by the rapid buckyball rotation. This motion also elim- 
inates any inter-Ceo J coupling 1 1 ] but leaves the dipolar cou- 
pling between spins on different buckyballs. Fig. [TJa) and 
[He) are the strongest evidence yet that pure dipolar coupling 



is sufficient to produce the phenomena 11211 . Moreover, dilu- 
tion of the spins is not required since ^^Y (/ = 1/2, 100% 
n.a.) experiments in Y2O3 [Fig.lTJb) andlHd)] also look quite 
similar to the earlier •^^Si (/ = 1/2, 4.67% n.a.) experiments 
in doped silicon ||5i, lZJ . 

The NMR signal in both Hahn echo and CPMG experi- 
ments is (x{Iyj,{t)) = X]f Tr{p(i)/j,;} in the rotating frame 
HI 0, B 01 • The time-dependent density matrix p{t) is cal- 
culated by starting with its conventional equilibrium value 
p(0) — = /zT' which assumes both the strong field 

and the high temperature approximations. Treating a strong 
90x pulse as a perfect ^ rotation about X, p{0) becomes 
p(0+) = e'(V2)/x^p(o)e-H'^/2)^x^ = /, 
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In between pulses, the spin Hamiltonian for these samples 
has two main parts. The first is the secular part of the dipolar 
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FIG. 1 : (color online). NMR measurements of ^"'C in Ceo (a), (c) and 
*^Y in Y2O3 (b), (d). Experimental parameters: T^SOO K, B=12 
Tesla, spin-lattice relaxation time Ti=26 s (Ceo), 2.3 hr (Y2O3). The 
dots in (a) and (b) are Hahn echo peaks while the red lines are CPMG 
echo trains. Simulated dipolar decay curves are shown for compari- 
son (black and blue, see text). Note the long tail for short interpulse 
spacings (a) and (b) and the even-odd effect for long interpulse spac- 
ings (c), (d). 
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where Bi 



[1 — 3cos^6'y ], 7 is the gyromagnetic 



the vector between spins i and j, satisfies 
Tij COS 6*^ (a static lab field points along z). The 

N 



ratio, and r, 

rij ■ z 

second part is a Zeeman term, Hn,, —^"—^zilzi, where 
Wz- is the angular frequency offset for spin i relative to on- 
resonance spins We further simplify this part 

by dropping the index i, since the linewidth of spectra stud- 
ied here suggest that cj^i is extremely uniform and entirely 
due to bulk diamagnetism of the powder [e.g., at B=12 Tesla, 
the linewidth r= 260 Hz (Ceo), 180 Hz (Si:Sb), and 3.1 kHz 
(Y2O3)]. The spin Hamiltonian describing free evolution is 
thus 



'Ho = Hzz + "Hn^ = Hzz + ^zlz 



(2) 



Two different approximations are used to calculate the ex- 
pected echo decays shown in Fig. [T] Both average over many 
disorder realizations (DRs), where each DR uses a randomly 
oriented lattice, with sites randomly occupied to match the 
natural abundance. In the first approximation, we build a 
closed quantum spin system around a central spin and keep 
all the terms in Eq. ([T]i while setting Viz — 0. Calculation 
of the evolution of the central spin is exact using numerical 
diagonalization [blue curve. Fig. [Tfa), with = 9 and 18 
DRs]. Unfortunately, computer limitations make it impracti- 
cal to include > 9 spins in this approach, so it fails for 
large or dense spin systems. In the second approximation, 
we drop the flip-flop terms in Eq. ([T]l, which yields an an- 
alytic expression ||5|] for any N (black curves. Fig. [T] with 
N 1000 and ^ 150 DRs). This second approach is unjus- 
tified for such clean samples, but Y2O3 is beyond the limits 
of the first method. Figure[T]shows Hahn echo data consistent 
with the expected decay due to the dipolar coupling. However, 
the CPMG echoes are observed well beyond this limit. 

The conventional (5-function pulse approximation [111 |2> HI 
0] treats very strong pulses as instantaneous vr rotations. In 
this limit, the density matrix for the Hahn echo (SEl) must 
agree with that for the nth Spin Echo (SEn) of a CPMG ex- 
periment. Thus the difference between the CPMG and Hahn 
echo data (Fig.[T]l is surprising jQl- 

Moreover, the experiments defined in Table |T] should pro- 
duce identical if the tt pulses are instantaneous. 
However, the measured echo trains for short r [Figs. |2fa) 
and|2jc)] exhibit a striking pulse sequence sensitivity (PSS). 
In this limit, the measured signal can either stay very close 
to one (CPMG, APCP), or dive rapidly towai'ds zero (CP, 
APCPMG). The PSS is observed UJ even when BBl com- 
posite TT pulses II14II are used. 

It is natural to attempt to blame the discrepancies between 
theory and experiment on extrinsic imperfections of instanta- 



TABLE I: tt pulse and echo phases for four pulse sequences of the 
form: 90x-{r-180^i -r-SEl-r-18O02-r-SE2-}''°P'=^*. 
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neous 180y pulses. Examples include misadjustment of the 
rotation angle, rf inhomogeneities, and pulse phase transients 
m [l]. We investigated these and other extrinsic errors of 
instantaneous pulses, but despite experimental improvements 
lfl3il . the effect persists. This led us to consider the limits of 
the conventional 5-function pulse approximation. 

The full Hamiltonian during an ideal pulse along the di- 
rection is 



Hp.— —hwilch 



-Ho 



(3) 



where the pulse's angular frequency, oji, is the same for all 
spins. The ^-function pulse approximation produces an in- 
stantaneous TT rotation about axis <j)i by keeping ujitp = tt 
while wi ^ 00 and tp 0, where tp is the pulse duration. 

In any real experiment, tp > 0, so the Hq term in Eq. (O 
might have some non-negligible effect. This possibility was 
first raised by Dobrovitski's simulations liisll of our silicon 
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FIG. 2: (color online), (a), (c) ^^Si measurements in Si:Sb (lO" 
Sb/cm^) with the four phase choices from Table|I] (b), (d) N — 7 
simulations for the same experimental conditions: ip=14 ^s, t—36 
Ais, r = 300 K, and 5 = 11.75 Tesla {^=99.5 MHz, ^ = 35.7 
kHz). Each curve is averaged over 1000 DRs with Qz/h drawn from 
a 290 Hz wide gaussian and typical couplings \Bi2\/h ~ 44.5 Hz, 
|_Bi7|//i « 3.5 Hz. Inset (a) shows rapid decay of CP data. Inset (b) 
shows a distinction between CP and CPMG simulations. The black 
dashed reference (a)-(d) sets = and turns off Ti.zz during each 
pulse. 
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data (3\. This is an intrinsic deviation from the instantaneous 
pulse Hmit, which cannot be avoided for tp > 0, as long as 
TCo is nonzero during the pulse. In fact, this is the sole cause 
of the small PSS that is seen in the N — 7 exact calculations 
[Fig Etb) and|2d)] since we have not included any extrinsic 
pulse errors. The dashed reference shown in Fig. |2ta)-|2|d) is 
the A^ = 7 exact calculation of CPMG for il^ = 0, with Ha 
set to zero during the pulses |16]. Given the strength of the 
pulses, — 108, we were surprised to see any difference in 
the simulations. 

To better understand the origin of the PSS in the simula- 
tions, we applied average Hamiltonian theory to CPMG with 
nonzero pulse duration (the other sequences will be treated 
elsewhere 11131 llal l. The toggling frame Hamiltonian 7i{ti) 
in Table In] is used to calculate the leading terms of the aver- 
age Hamiltonian {H = 7i(") + U'-^^ + . . .) fl B H- The 
zeroth order term is 
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CPMG ^ -t (4'''^z 



(4) 



where the cycle time tc — 4r+2ip. The time-dependent terms 
within each pulse (see Table HHi give rise to a first order term: 
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In contrast, the full H for CPMG with (5-function vr pulses is 
simply Hzz- Since the experimental consequences of finite 
pulses are not obvious from Eqs. (|4|i and (|5]l, we used simu- 
lations to study their effects. Focusing on the simplest case of 
= leaves only one commutator in Eq. (|5]l. 
Figure [3] shows calculations of the CPMG sequence that 
yield a long tail with parameters ilz—0, = 40 kHz, and 
T = 1 /is. Setting Ha to zero during the tt pulses yields a 
leading term ^Hzz, causing the fastest signal decay [Fig. 
Qla.)]. Keeping Ha during each pulse modifies this leading 
term to Hqpj^q slowing the signal decay [Fig.[3|b)]. Adding 
in the higher order correction (Hqpmq) slows the decay even 



TABLE II: Toggling frame Hamiltonians Ti.{ti) during each inter- 
val i for the CPMG cycle {r-180y-2r-180y-r} with pulse time tp. 



Here, H,,=Ef>, B,, (s/, J,, -/.•/,), W^^=Ef>, fS,, (1.,/.^ + 
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hJ^^), and = Ef>, iBrjihJ,^ - The factors 

(Ce-, Se^, C26i, S2ei)=(cos9i, sin 9i, cos20i, sin 2^0 are time- 
dependent, as = Luiti, and < ti < Ti |3, 4, ^. 
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more, resulting in the long tail [Fig. [3fc)]. This is true even 
for^z + Hi. 

Average Hamiltonian theory was first used to design line- 
narrowing sequences by generating particular leading terms in 
the average Hamiltonian. For example, pulse sequences can 
be designed to set H^'^^ =0, which results in a lack of signal 
decay. Higher-order corrections would then modify 

this behavior by causing some decay in the signal aaa. 
In contrast. Fig. [3] shows that for the CPMG sequence, higher 
order corrections to the zeroth order Average Hamiltonian can 
slow decay and even cause a long tail to develop. 

The exact calculations in Fig. |3ld)-[3f) look very similar at 
early times, as expected for the conditions of the simulations 
[Fig. [3linset)]. Most surprisingly, the tail height grows with 
^ [Fig- [3d)-[3if)]. Thus, knowing the linewidth or By scale 
is not enough to predict the shaM of the whole curve. Extrap- 
olating from these results |I3lll6ll . it appears that the small 
and Bij used in Fig. |2b) inhibit any tail. Moreover, a notice- 
able tail should emerge in Fig. |2|b) if the simulations could 
use N > 20. 

While the long-lived coherence in CPMG is reminiscent of 
spin locking [3], the underlying dynamics caused by strong tt 
pulses induce important time-dependent changes in the sys- 
tem. More information is revealed by visualizing the entire 
time-evolved density matrix as shown in Fig. H) After two 
strong TT pulses (Fig. 01 leftmost column) the density matrix 



still looks like the initial state p(0'^) 



The top two 



§ 60- 

-e 40 

20 

0- 





mi/2it.40kHz 


I. ■ 







-20 -10 10 20 30 40 
Relative Frequency (kHz) 




Force //„ =0 
during ;i -pulses 



(b) N=4 
' "(a)" N-4" " 



10 15 
Time (msec) 



FIG. 3: (color online). CPMG calculations for pure dipolar decay 
(fiz = 0, =40 kHz, T— 1/is). Each curve averages 400 DRs 
[exception (f): 80 DRs] of N spins on a silicon lattice, with ^'=5"/ of 
^''Si. Several approximations are used for A'^=4: (a) setting Ho ~ 



during pulses, (b) using H 



(0) 



only, and (c) using H 



H 



(1) 

CPMG 



(0) 

CPMG 



+ 

only. Exact calculations for (d) A'' = 4, (e) N — 6, and (f) 
N=8 show that the tail height depends on A'^ (even N are compared 
to avoid artifacts figll ). Inset shows the pulse strength used (lui 3> 
2ttT) compared to the calculated spectra for A =4 (red) and A =6 
(blue). 
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FIG. 4: (color online). Calculated snapshots of a 6 spin (2® x 2^) den- 
sity matrix evolving during CPMG with conditions as in Fig. |3] The 
red-white-blue color scale shows the phase angle; black cells have 
negligible magnitude. The top two rows use the 5-function pulse 
approximation. The bottom two rows use exact tt pulses. t2Q1 



rows use the (5-function pulse approximation. In each DR (first 
row), the flow of coherence is restricted to the ±l-quantum 
coherence cells |ll'4l], consistent with the selection rules for 
pure dipolar coupling and ^-function tt pulses. Over 400 DRs 
(second row), all quantum coherences average to zero after 
many pulses, so the measurable signal decays as in Fig. 13 a). 
The bottom two rows use exact tt pulses [see Eq. ^]. In each 
DR (third row), coherence flows throughout the entire den- 
sity matrix, since exact tt pulses [see Eqs. (|4|i and (|5]i] open 
many coherence transfer pathways [4]. Despite this complex- 
ity, observable coherences emerge after averaging over 400 
DRs (fourth row), even after many pulses, a counterintuitive 
result that yields the smooth curve in Fig.[2e). 

Nonzero pulse duration has been studied for hne-narrowing 
sequences built around 7r/2 pulses, and it is not considered to 
be an error For example, the leading term of the average 
Hamiltonian for the Ostroff-Waugh experiment iHtIi is iden- 
tical for both instantaneous and nonzero duration 7r/2 pulses. 
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In contrast, nonzero duration tt pulses 



introduce completely new operators (^ Hzz) into the aver- 
age Hamiltonian of CPMG [Eqs. (HJl and ©J. Nonzero du- 
ration effects should be maximal for m x tt pulses, where 
m — 1,2,3,.... 

In general, related effects could arise whenever the applied 
pulse term does not commute with the spin-spin interaction 
(e.g., Ising or general anisotropic Heisenberg couplings). Ex- 
amples include ESR experiments on dilute moments or bang- 
1 — I 

bang control sequences [18] applied to systems of qubits with 
weak always-on coupling. Local pulses that address only a 
subset of coupled spins are also susceptible to these effects. 



Developing an improved understanding of these many-body 
corrections to pulse action will enable the rational design of 
pulse sequences optimized to boost signal-to-noise, narrow 
spectra, and achieve desired coherence transfer pathways. 
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